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1. Introduction 



Accurate perturbative calculations beyond leading order in quantum chromo dynam- 
ics (QCD) are an important ingredient in improving our understanding of jet pro- 
duction in current and future high energy collider experiments at the Tevatron and 
LHC. At present, next-to-leading order calculations have become standard and are 
used to make comparisons with experimental data. For example, the next-to-leading 
order O (a^) predictions for jet production in pp collisions [1, 2] based on the one-loop 
matrix elements computed by Ellis and Sexton [3] have been successfully compared 
with a wide variety of experimental observables using data from the Tevatron and 
the CERN SppS. To date these comparisons have been limited by both experimental 
and theoretical uncertainties at the 10% level. However, improvements in detector 
technology, as well as the expected large increases in the luminosity of the colliding 
particles, should significantly improve the quality of the experimental data and will 
require more accurate theoretical calculations either to claim new physics or to refine 
our understanding of QCD. 

The theoretical prediction may be improved by including the ncxt-to-ncxt-to- 
leading order perturbative predictions. This has the effect of (a) reducing the rcnor- 
malisation scale dependence and (b) improving the matching of the parton level 
theoretical jet algorithm with the hadron level experimental jet algorithm, because 
the jet structure can be modeled by the presence of a third parton. The improve- 
ment in accuracy expected at next-to-next-to-leading order can be estimated using 
the renormalisation group equations together with the existing leading and next-to- 
leading order calculations and is at the 1-2% level for centrally produced jets with a 
transverse energy. Et. of around 100 GeV. 

The full next-to-ncxt-to-leading order prediction requires the knowledge of the 
two-loop 2-^2 matrix elements as well as the contributions from the one-loop 
2-^3 and tree-level 2^4 processes. At large transverse energies, Et ^ ^^iquark, 
the quark masses may be safely neglected and we therefore focus on the scattering 
of massless partons. Techniques for computing multiparticle tree amplitudes for 
2 — > 4 processes, and the associated crossed processes, are well understood. For 
example, the helicity amplitudes for the six gluon gg — > gggg, four gluon-two quark 
qq — > gggg, two gluon-four quark qq q'q'gg and six quark qq q'q'q"q" have been 
computed in Refs. [4, 5, 6, 7]. Similarly, amplitudes for the one-loop 2 — 3 parton 
sub-processes gg — >■ ggg, qq — > ggg, qq — > q'q'g, and processes related to these by 
crossing symmetry, are also known and are available in [8, 9, 10] respectively. 

Although the two-loop contribution for gluon-gluon scattering in = 4 super- 
symmetric models has been known for some time [11], the evaluation of the two- loop 
2^2 contributions for QCD processes has been a challenge for the past few years. 
This was mainly due to a lack of knowledge about planar and crossed double box 
integrals that arise at two-loops. In the massless parton limit and in dimensional 
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regularisation, analytic expressions for these basic scalar integrals have now been 
provided by Smirnov [12] and Tausk [13] as series in e = (4 — -D)/2, where D is 
the space-time dimension, together with algorithms for reducing tensor integral to a 
basis set of known scalar (master) integrals [14, 15]. This makes the calculation of 
the two-loop amplitudes for 2 — > 2 QCD scattering processes possible. 

Following on from the pioneering work of Bern, Dixon and Ghinculov [16] who 
completed the two-loop calculation of physical 2 — > 2 scattering amphtudes for the 
QED processes e^e" — > //"^//~ and e"^e~ — > e~e'^ , we have studied the O (a^) contri- 
butions arising from the interference of two-loop and tree-level graphs for the QCD 
processes of quark-quark [17, 18, 19] and quark-gluon [20] scattering. In these pa- 
pers we presented analytic expressions for the infrared pole structure (that ultimately 
cancels against contributions from the 2 — 3 and 2 — > 4 processes), which agrees 
with that anticipated by Catani [21], as well as the finite remainder. 

To complete the set of two-loop contributions to parton-parton scattering re- 
quires the study of (non-supersymmetric) gluon-gluon scattering. Bern, Dixon and 
Kosower [22] were the first to address this process and provided analytic expressions 
for the maximal-helicity- violating two-loop amplitude. Unfortunately, this amplitude 
vanishes at tree level and docs not contribute to 2 — >^ 2 scattering at next-to-next- 
to- leading order 0{af). It is therefore the goal of this paper to provide analytic 
expressions for the O (af ) two-loop corrections to gluon-gluon scattering 

g + g-^g + g. (1.1) 

As is in Refs. [17, 18, 19, 20], we use the MS renormahsation scheme to remove 
the ultraviolet singularities and conventional dimensional regularisation. where all 
external particles are treated in D dimensions. We provide expressions for the in- 
terference of tree-level and two-loop graphs. The infrared-pole structure agrees with 
that obtained using Catani's general factorisation formulae [21]. The finite remain- 
ders arc the main new results presented in this paper and we give explicit analytic 
expressions valid for the gluon-gluon scattering process in terms of logarithms and 
poly logarithms that are real in the physical domain. For simplicity, we decompose 
our results according to the powers of the number of colours N and the number of 
light-quark flavours Np- 

Our paper is organised as follows. We first establish our notation in Sec. 2. 
Analytic expressions for the interference of the two-loop and tree-level amplitudes 
are given in Sec. 3. In Sec. 3.1 we adopt the notation used in Ref. [21], to isolate 
the infrared singularity structure of the two-loop amplitudes in the MS scheme in 
terms of the one-loop bubble integral in D = 4 — 2e and the one-loop box integral in 
D — 6 — 2e. Analytic formulae connecting these integrals in the various kinematic 
regions are given in Appendix A. We demonstrate that the anticipated singularity 
structure agrees with our exphcit calculation. The finite C (e°) is given in Sec. 3.2 
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in terms of logarithms and polylogarithms that have no imaginary parts. Finally we 
conclude with a brief summary of the results in Sec. 4. 



2. Notation 

For calculational purposes, the process we consider is 

9{Pi) + 9{P2) + 9{P3) + 9{pa) 0, (2.1) 
where the gluons are all incoming with light-like momenta, satisfying 

pt+P^+P^+p't^O, Pi^O- 
The associated Mandelstam variables are given by 

S = (P1+P2)', t^{p2+P3f, u^{pi+p3f, s + t + u^O. (2.2) 

The gluons also carry colour indexes, Oj, in the adjoint representation. 

We work in conventional dimensional regularisation treating all external quark 
and gluon states in D dimensions and renormalise the ultraviolet divergences in the 
MS scheme. The renormalised four point amplitude in the MS scheme can be written 



\M) = 4:7ras 



(2.3) 

where ag = q;s(//^) is the running coupling at renormalisation scale // and the |A^^'^) 
represents the colour-space vector describing the renormalised i-loop amplitude. The 
dependence on both renormalisation scale /j, and renormalisation scheme is implicit. 
We denote the squared amplitude summed over spins and colours by 

{M\M) ^J2\M(9 + 9^9 + 9)\'^ V{s, t, u). (2.4) 

which is symmetric under the exchange ois,t and u. The function T> can be expanded 
perturbatively to yield 

.2% J ^"'"'"^ ' V27r. 



V{s,t,u) = IGn^al 



V\s, t, u) + (^) V%s, t,u)+(^y V'{s, t,u) + (at) 



(2.5) 

where 

V\s,t,u) = 

^16mi-e)^(3-^-^-i^), (2^6, 

V%s,t,u) = [{M^^^\M^'^) + , (2.7) 

V%s,t,u) = + (A^(o)|A^(2)) + (M(2)|M(o))) ^ (2.8) 
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where N is the number of colours and V = N"^ — 1. Expressions for T)^ are given 
in Ref. [3] using dimensional regularisation to isolate the infrared and ultraviolet 
singularities. 

In the following sections, we present expressions for the infrared singular and 
finite two-loop contributions to 

= + (2.9) 

We defer the self-interference of the one-loop amplitudes 

V^^^''^\s,t,u) = (2.10) 

to a later paper. 

As in Refs. [17, 18, 19, 20], we use QGRAF [23] to produce the two-loop Feynman 
diagrams to construct We then project by {M.^^^\ and perform the summa- 

tion over colours and spins. It should be noted that when summing over the gluon 
polarisations, we ensure that the polarisations states are transversal (i.e. physical) 
by the use of an axial gauge 

^ efe- = -g^^- + ^P^'' + ^^i' (2.11) 

spins P'^ 

where pi is the momentum, Cj is the polarisation vector and rii is an arbitrary light- 
like 4-vector for gluon i. For simplicity, we choose rii = P2, ni^ = Pi, n-^ = P4 and 
^4 = Pa- Finally, the trace over the Dirac matrices is carried out in D dimensions 
using conventional dimensional regularisation. It is then straightforward to identify 
the scalar and tensor integrals present and replace them with combinations of the 
basis set of master integrals using the tensor reduction of two-loop integrals described 
in [14, 15, 24], based on integration-by-parts [25] and Lorentz invariance [26] identi- 
ties. The final result is a combination of master integrals in D = 4 — 2e for which 
the expansions around e = are given in [12, 13, 14, 15, 24, 27, 28, 29, 30]. 

3. Two- loop contribution 

We further decompose the two-loop contributions as a sum of two terms 

p8(2xo)^^^^^^^ ='Poles{s,t,u) + Tinite{s,t,u). (3.1) 

Voles contains infrared singularities that will be analytically canceled by those oc- 
curring in radiative processes of the same order (ultraviolet divergences are removed 
by renormalisation) . ^Finite is the remainder which is finite as e — > 0. 
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3.1 Infrcired Pole Structure 

Following the procedure outlined in Ref. [21], we can write the infrared pole structure 
of the two loop contributions renormaliscd in the MS scheme in terms of the tree 
and unrenormalised one-loop amplitudes, |A^*^°)) and lA^^^'""^)) respectively, as 



Voles = 2 Re 



+ (M(°)|lW(e)|M(^'"")) 
+ (M(°)|i/(^)(6)|M(°)) 



(3.2) 



where the Euler constant 7 = 0.5772 The first coefficient of the QCD beta 

function, Pq, for Np (massless) quark flavours is 

^^^ nc.-mN. _ c.^iv, r«a, (3.3) 

and the constant K is 

Note that the unrenormahsed one- loop amplitude is what is obtained by 

direct Feynman diagram evaluation of the one-loop graphs. 

It is convenient to decompose and \A4^^''^"-^) in terms oi SU{N) matrices 

in the fundamental representation, T", so that the tree amplitude may be written 
as [31, 32, 33] 

]X(°)) = Tr(r»ir"2T«3T»4)^t/^^(l,2,3,4), (3.5) 

P(2,3,4) 

while the one-loop amphtude has the form [34, 35] 

l_y^(i,«n)^ ^ ^ j2 Tr(r"ir"2r"3r"^)>ii^l(i,2,3,4) 

P{2,3,4) 

+ Tr(T"iT"2)Tr(T"3T"4)y^i;!3(l,2,3,4) 

Q(2,3,4) 

+ Np J2 Tr(r»ir"2r«3r"*)y4i;f^(i,2,3,4). (3.6) 

P(2,3,4) 

In these expressions Zlp(2,3,4) ^^^^s over the 6 permutations of indices of gluons 2, 3 
and 4 while J2q(2,3,4) includes the three choices of pairs of indices, as it is further 
detailed in Eq. (3.9). We note that the tree subamphtudes are further related by 
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cyclic and reflection properties as well as by the dual Ward identity [32, 36] and 
more general identities [37], while the subleading-colour loop amplitudes ^4^3 are 
related to the leading-colour amplitudes ^4:'^ [34, 35]. Some of these relationships 
are made explicit using an alternative basis in terms of SU (N) matrices in the adjoint 
representation [38]. 

To evaluate Eq. (3.2) we find it convenient to express \A4^^^) and as 
nine-dimensional vectors in colour space 

= (Ti, Ta, %, Ts, %, 0, 0, of , (3.7) 
\^{i,un)^ = (A, £2, A, A, A, A, Cr, Cs, C,f, (3.8) 

where Q'^ indicates the transpose vector. Here the 7^ and £j are the components of 
and 1^(1'^^")) in the colour space spanned by the (non-orthogonal) basis 

£7^ _ rYj. ^rpairpa2rpa:irpa4-j 
_ 'Jij, ^rpairpa2rpa4rpa3^ 

_ ^J^airpa3rpa2rpa4-\ 

rj^^ I^rpairpa3rpa4rpa2^ 

_ ^rpairpa4rpa:irpa2^ 

C7 = Tr (r"iT"2) Tr (T^^T"^) , 
Cg = Tr [T^^T"^) Tr [T^^T"^) , 

Cg = Tr (T"ir"^) Tr (T^^T"-^) . (3.9) 



The tree and loop amplitudes % and Ci are directly obtained in terms of A. 



tree 



^4^4, ^4^3 and .4.4^1^^ by reading off from Eqs. (3.5) and (3.6). As we will see, the 
amplitudes themselves are not required since we compute the interference of tree and 
loop amplitudes directly. 

In the same colour basis, the infrared-singularity operator l'^^\e) introduced by 
Catani [21] has the form 



l[l/2] 



X 



( N{S 







(S-U) 


V (T-U) 






N{S + U) 
Ar(T + U) 








(S-T) 
(U-T) 










(U-S) 
(T-S) 







iV(T + U) 






(u 

(T 



s) 
s) 









iV(S+U) 


(S-T) 
(U-T) 












(T-U) (S - U) \ 
(U - T) (S - T) 
(T - S) (U - S) 





(U-T) 



iV(S + T) (T - U) 
(S - U) 2iVS 


(T-U) 



(T-S) 
(S-T) 





2A^U 




(U-S) 


(S-U) 



2ArT / 
(3.10) 
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where 



/f2 

s 



' t 



U 



u 



(3.11) 



The matrix I^'^^ (e) acts directly as a rotation matrix on and in colour 

space, to give a new colour vector \X), equal to I^^\e)\M^^^), I^^\e)I^^\e)\M^^^) or 
jW(e)|Al(i'™)). 

The contraction of the colour vector \X) with the conjugate tree amplitude obeys 



the rule 



{m^'^\x) = y: e e vxjc:c, 

spins colours i,j=i 



(3.12) 



In evaluating these contractions, we typically encounter Z^coiours Q which is given 
by the ij component of the symmetric matrix OC 

( Ci C2 C2 C2 C2 C3 -A^ \ 

C2 Ci C2 C2 C3 C2 NV NV -N 

C2 C2 Ci C3 C2 C2 -N NV NV 

C2 C2 Cs Ci C2 C2 -N NV NV 

C2 C3 C2 C2 Ci C2 NV NV -N , (3.13) 

C3 C2 C2 C2 C2 Ci NV -N NV 

NV NV -N -N NV NV N'^V N"^ N"^ 

-N NV NV NV NV -N N'^ N'^V N^ 

\ NV -N NV NV -N NV N'^ N"^ N'^V j 



CC 



V 



16A^2 



with 



Ci = N^- 3Ar2 ^ 3^ 



3 - TV^^ Ca = 3 + TV^ 



(3.14) 



Similarly, we find that the interference of the tree-level amplitudes X^, 
given by TTj^ , where 



T*T is 

spins -'i -^3 



TT 



QA{1- ef{t^ + ut + u 



2\2 



and the vector V is 



V 



(m, t, s, s, t, M, 0, 0, 0) , 



(3.15) 
(3.16) 

while the interference of the tree-level amplitudes with one-loop amplitudes Z^gpins ^C^j 
is given by TTy, where 

rC = V^W, (3.17) 

and the vector W is 

w^[T{s,t), T{s,u), j^(u,t), j^(u,t), j^(s,u), j^(s,t), g, g, g). (3.18) 
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Here the function JF(s, t) is symmetric under the exchange of s and t, while Q is 
symmetric under the exchange of any two Mandelstam invariants, so that 

JF(s, t) = his, t, u) + hit, s, u), (3.19) 
G = his, t, u) + his, u, t) + hit, s, u) + hit, u, s) + hiu, s, t) + hiu, t, s). 

(3.20) 

Here fi and /2 are given in terms of the one-loop box integral in D — 6— 2e dimensions 
and the one-loop bubble graph in D = 4 — 2e, 



his,t,u) 



+ 



+ 



IGNil - 2e) 

^¥ 

SNpil - 2e) 

7t 

l6Ar(i -6) 

52^2^^6(3 - 2e) 

+ (36 - 99e + 93e^ - 24e^ - 2e^) sH'^ + (1 - e) (24 - 50e + 23e^) st^ 

+4(1 - e)(l - 2e)(3 - 2e)t^] Bub(t) 

16Nf 
st^uiS - 2e) 



2(1 - {s^ + s^t + st^ + t^) + 3(1 - 5e)sH^] Box%s, t) 
il-ef (s^ + i^) + e(l + 3e)si] Box^(s, t) 
(l2 - 22e + 12€^ + 2e^) / + (24 - 58e + SOe^ - - 2e^) sH 



(4 - 12e + 166^ - 4e^) + (3 - lOe + 23e^ - 8e^) sH 



+ (6 - 15e + 21e^ - Se'M + (1 - e) (5 - 6e + 2e^) Bub(t), (3.21) 



, 32(1 -2e) 
his,t,u) = ^ , ' 

32(1 -e) 



esu^ 



-4(1 - e)^si + 3(1 - 5e)ii2] Box^(ii, i) 
4(1 - 2e)(l - e)t2 + (8 - 17e)(l - e)ut 



+ (6 - 20e + 15e^ + e^) u'^] Bub(s). 



(3.22) 



Series expansions around e = for the one-loop integrals are given in Appendix A. 

Finally, the last term of Eq. (3.2) that involves H^'^\e) produces only a single 
pole in e and is given by 



3^7 



(3.23) 



where the constant i/^^^ is 

(2(3 + - + - 

V ^ 3 36 



20 



TT^ 89 



tt" AT^ + — iV^ + — - - — NNp - -z^ 



27 



18 27, 



A^ 



(3.24) 



and Cn is the Riemann Zeta function with (2 = 7r^/6 and Cs = 1.202056 . . . We note 
that if(2^ is renormahsation-scheme dependent and Eq. (3.24) is valid in the MS 
scheme. We also note that Eq. (3.24) differs from the corresponding expressions found 
in the singularity structure of two-loop quark-quark and quark-gluon scattering. This 



8 



is due to double emissions from the gluons. In fact, we note that H'^'^^ for quark- 
gluon scattering is the average of the H^"^^ for gluon-gluon scattering and quark-quark 
scattering, as may be expected by counting the number of different types of radiating 
partons. 

It can be easily noted that the leading infrared singularity in Eq. (3.2) is O (1/e''). 
It is a very stringent check on the reliability of our calculation that the pole structure 
obtained by computing the Feynman diagrams directly and introducing series expan- 
sions in e for the scalar master integrals agrees with Eq. (3.2) through to 0{l/e). 
We therefore construct the finite remainder by subtracting Eq. (3.2) from the full 
result. 

3.2 Finite contributions 

The finite two-loop contribution to V^{s,t,u) is defined as 

J^inite{s, t,u) ^ V^^'^''^\s,t,u) - Voles{s, t,u), (3.25) 

where we subtract the series expansions of both (2x0) ^) ^j^j Voles{s, t, u) and 
set e — s> 0. As usual, the polylogarithms Li„(-u;) are defined by 



Li„(w) = / — Li„_i(t) forn = 2,3,4 

JO t 

fW rjf 

U,{w)^- ^log(l-i). (3.26) 
Jo t 



Using the standard polylogarithm identities [39], we retain the polylogarithms with 
arguments x, 1 — x and {x — l)/x, where 

t 1J- u X — 1 
x = — , — = l-x, z^--^ . (3.27) 

S S t X 

For convenience, we also introduce the following logarithms 

X = log(Z^), r = log(^), 5 = log(-i), (3.28) 

where /j, is the renormalisation scale. 

We choose to present our results by grouping terms according to the power of 
the number of colours N and the number of light quarks Np, so that 

Finite(s, t,u)^V (^N^A + N^B + N^NpC + NNpD + N^Nj,E + N^F^ , (3.29) 

where 

^ = I U8 Li4 (x) - 48 Li4 (y) - 128 Li4 (^) + 40 Lig (x) X - 64 Lig {x)Y Lig {x) 
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98 16 

+64 Li3(y) X - 40 Li3(y) Y + 18 Us(y) + - Li^ix) X-- Li^ix) ^ - 18 U^iy) Y 

_^X^ + 28X3F--X3-16X2F2 + l^X2y_ 35^2^2_ 38^2_ 22^^2 
6 3 3 3 3 3 

20 31 22 37 

— xy=^-9xy2 + 8xy7r^ + ioxy — Xt^"^ -^nux^ — sx^—x 

3 12 3 27 

11 ,,4 41 . 11 ^2 2 22 ^^. 266 ^. 35^^ 2 418^^^ 257^^ 

H — — Y^ — y^TT^ — SY^ ^ y^ — yvr^n sy a y 

6 9 3 3 9 12 9 9 

^18/^ V 31 4 11c, 2^31 2^ 242 2 ^418 2156 

+ 18 Csy TT O TT H TT H 6 H H -J 

30 9 9 9 9 ^ 27 

/ 64 
+ ( - 256 Li4(x) - 96 l^^iy) + 96 Li4(^) + 80 l^\z{x) X + 48 Li3(x) - y Lisl^) 

-48 Li3(2/) X + 96 Li3(|/) y - ^ Li3(|/) + y U^{x) X - | Li2(a;) vr' + ^ Lisd/) Y 

+^ _ ^ x3 y - ^ x3 + 20 y 2 + ^ y + 24 X + 76 - ^ 5 
+^xy3 + l^xy^-i^xy.^ + iZ^^xy-i^xy-^x.^-48C3X 

3 3 3 3 3 3 

2350 ^ 440 ^^ 176,^3 4 ^ „ 2 176^,^2 494^^ 2 5392 ^^ 

+^x + — 5x + 4y^-— y3 + -y2 7r'-— ,5y'-— yvr' + ^y 

, , 496 4 308 ^ 2 200 2 968 ^2 8624 „ 44372 
9 j u 

+ Li3(x) - ^ Li2(x) X + 2 x" - sx" y - ™ + 12 X' + ^ x= r + 5 X' 

_|sx^ + ^x^_8xr3_^xy.^ + lp5xy-fx.^ + H5x 

- 8C3X + 4y' - H5y'- f yV= - H5sy= _ !?5y,^ _ 16C3y 

5392 ^^ 4 4 308 ^ 2 ^^2 ^ . 1408 ^ 968 ^2 44372 
8624 \ ^2 

+ f- Li3(x) - - Li2(x) X - + — X3 - - ^2 y + - ^2 TT^ + - 5 

\^3 ^^3 ^^ 9 3 3 3 

_l^x^_ioxy + ilA-.= + 4c,x-l?lsx-™x + l.'-51.= 

9 2 ^ 9 27 30 9 

11 2 418^ 242 . 2156^ 11093\Mt 
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+ 1^-176 Li4(x) + 88 Li3(j;) X - 168 Li3(x) ^ - ^ Li3(x) + ^ Li2(x) X 

''^ SX-' I^AxyJ+?^xy + usxl -mx\^ + ?^l 



3 3 3 36 27 

1870 ^,, .o. 31 4 361 2 517 ^ ^ 1331^2 12452 ^ 

+Hi!c3-^-44SC3)UU4. (3.30) 



I - 288 LU{x) + 480 Li4(?/) - 288 Lu{z) + 240 Li3(x) X - 144 Li3(x) F 

+224 Li3(a;) + 144 U^iy) X - 432 U^iy) Y - 224 U^iy) + 48 Li2(x) 

-224 Li2 (x) X - 176 Lis (a;) tt^ + 48 Lis (y) + 224 Lis (y) F - 16 + 1 12 F 

X^ - 48 X^ r 2 + 180 X^ r - 40 X^ TT^ + 220 X^ - 32 X _ 92 X 

3 

376 9Q9 
-16XY7r^-^XY-16X7r^-S0X + 96CsX + 8Y^ + —Y^-32Y^7r^ 
3 3 

y 2 ^ 16 y TT^ + 80 y - 96 C3 1^ + - TT^ - 18 TT^^ 



+ ^ - 576 Li4(x) + 384 Li4(y) - 1152 Li4(z) + 1056 Li3(x) X - 768 Li3(x) Y 

+448Li3(x) + 768Li3(y)X- 768Li3(y)y + 896Li3(y) - 192 Lis (x) 

-448 Lis (x) X - 544 Li2 (x) TT^ - 384 Li2 (y ) X y - 896 Lis (y ) r - 28 X^ + 144 X^ y 
320 

+ X^ - 336 X^ y 2 - 224 X^ y - 40 X^ tt^ - 64 X^ - 32 X y ^ + 128 X Y^ 

3 

1 888 

-64 X y + X y - 288 X + 160 X - 1248 (sX - 240 Y^ - 928 Y 

3 

. 1216 4 1912 2 ..o . \ ^ 
+768C3l^ + -r^7r^-^7r2-448C3 - 
15 3 J u 

+ 1^-384 Li4 (y) - 384 Li4 (z) + 384 Li3 (x) X - 384 Li3 (x) y + 384 Li3 (y) X 

-192 Li2(x) X2 - 192 Li2(x) - 384 Lis(|/) X y - 8 X^ - 32 X=^ y - 176 X=^ 

752 

-192 X2 y2 + 352 X2 y - 80 X2 + X^ - 32 X y tt^ - 176 X tt^ - 384 (3 X 

3 



-96 y 2 7r2 - 352 y 772 + 384 (3 1^ + 56 tt^ - — 772) ^ 

3 



+ 1^-192 Li4(x) + 192 Li3(x) X - 96 Li2(x) X^ - 4 X^ - 32 X'^ y + 88 X 
+12X2y2 -88X2 y + 48X2 TT^ _^X2-48Xy7r2 + ^Xy 
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64 4 -, o 2 \ ""^ 



+ 1^48 Li3(x) X + 144 Li3(x) Y + 672 Li3(x) - 48 U2{x) X' - 672 U2{x) X + IQX 
4 1444 

-32 r - - + 24 + 12 r - 192 tt^ + — - + 72 x r tt^ 

o o 

+— X y - 624 X TT^ + 80 X - 288 Cs X + — TT^ - 707 TT^ - 36 - 2800 Cs 
3 15 



)} 



+<^w^n, (3.31) 



C = 1 1| - 24 Li4(x) + 24 hu{y) + 112 hu{z) - 44 Li3(x) X + 56 Li3(x) ^ + y Li3(a;) 

74 32 

-56 Li3(y) X + 44 Li3(y) Y -22 U^{y) - - U^ix) X + - U^ix) tt^ + 22 U^iy) Y 

26X^^ + 4X3 + 14X2^2 F + 7X2 TT^ + —X^ + 5 5X2 

4 3 2 

22 , , , 31 , 637 26 

+— xy^ + iixy2 _4xy7r^ -iixy + — X7r^ + i2C3X x sx 

3 6 ^ 27 3 

12 9 3 18 3 6 9 

98 1 4 2 ^ 2 203 2 4 ^ 88 „ 4849 386 \ ^ 

/ 124 
+ ( 224 hu{x) + 48 hu{y) - 48 hu{z) - 88 U^{x) X - 24 Us{x) ^ + -3- Li3(x) 

280 124 64 

+24 Li3(y) X - 48 Li3(y) Y + — Us(y) - — U2(x) X + - Li^ix) 

_'^U2(y)Y--X^ + 6X^Y--X^-3X'Y'--X'Y--X^7r'-2SX' 
3 6 3 3 3 

70 , o , 2 „ 148 22 , 

-y X2 -6Xy^ -26Xy2 __XYTr^ + 4SXY+—XY-—XTr^ 

A V 829 , 10 J 356 2 352 , 1544 388 9698\ i 

^^^'-'^--^^ -T^" -17" ^ '^'+"8rj« 

\^ 3 ^ ^ 3 ^ ^ 9 3 9 3 3 

14^ 2 224 ^ 352 ^^ 32 „ 32 ^„ 116^^ 2 1096 ^^ 56^ 2 

340 2 1544 ^ 9698 32 ^ 352 ^2\ 

"^^+-8r + ¥^^"-9-^ h? 
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+ (^-^U,ix) + lu,ix)X-'fx' + ^X'Y-^-^X'-lsx' + nXY-Xn' 
112 „ 176 „ „ 203 o 88 4849 386 ^ 4: \ut 



+ (136 Li4(a;) - 68 Us{x) X + 120 Li3(a;) Y + ^ Us{x) - ^ U2{x) X-^X^ 



j,}^ X^Y -— X^ + lhX^Y^ + - X^Y ~— X^ + — X^ + — S X"" 
3 9 3 3 18 3 

3 6 18 27 9 

47 4 3799 2 47 ^ 2 484 2825 ^ 932 ^ 70025 \ 1 

+ (3-32) 



f / 46 
D = l\2A hu{x) - 24 hu{y) + 88 Li4(^) - 52 U^{x) X + 36 Li3(x) ^ - y Li3(x) 

46 46 44 

-36 Li3(2/) X + 52 Li3(2/) F + y Li3(y) - 4 Li2(x) + y Li2(x) X + — U^ix) tt^ 

-16 Li2(y) X y + 4 Li2(y) - ^ ^^(y) F + ^ - | y + ^ + 3 

184^2.. 13^2 2 545^2 38 136^^. . 155^^^ 

y H — H — X y^ H X y^ + - X y H x y 



3 3 6 3 3 3 3 

-iox7r2-32C3X + -x--y4-— y3 + 3y2 7r2 + — y2 + ioy7r2--y 

3 12 18 6 3 



2 



/ 92 
+ ( 176 Ui{x) - 48 hU{y) + 48 hu{z) - 104 Li3(x) X + 32 U^{x) Y 

-32 UM ^ + 64 Li3(?/) >^ - ^ Li3(l/) - 8 Li2(a;) + y Li2(a;) X + ^ U^ix) 

+16 Li2(y) X y - 16 Li2 (y) y' + ^ Li2(y) y - ^ - 10 X^ y - ^ X^ + 19 X^ y^ 

o 6 9 

+l^X2y-17X2 7r2 + -X2--Xy3-87Xy2--Xy7r2- 260Xy 
3 3 3 3 

+^X7r^-l^X + 168C3X + 7y^7r^ + ^y7r^ + 8y-32C3y-^7r^ 

742 2 188 ^ 110 „ \ t 
3 3 3 ) u 

+ f32X3-64X2y-^X2 + 32X7r2 + 64y7r2 + 8y + — 7r2 + 85 

y 3 3 
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+ 1^64 Li4(a;) - 20 U^{x) X - 108 U^{x) F - 46 Li3(a;) - 12 Li2(a;) 

+A%U,{x)X + ^X'-lQX^Y-^X^-'^-^X^Y^-^4x'Y-\^'^^ 
S.Z lo z o o 

1 _ 16 ^ y _ U y ^ 340 ^ ^ ^ ^ _ 52 ^ _ 67 



6 3 6 3 3 20 

981 2 o 1166 , 461 

TT^ + 11 5^ H Cs 

12 3 ^ 12 



2981 2 o 1166 , 461\ 1 f 1 
+-77^ TT^ + 11 + ^— C3--rz- \ + {u^t\, (3.33) 



1^33 33 3 3393 

+ -Y t:^ + - SY + - Y ^ 7r^ + -5^2 — 

3 9 3 27 9 y s2 

/2o2, 4,4,2 ,8 2,8 
+ { - X^ - - X^Y + - X^ + - S X^ - - XY^ - - S XY + - X Ti^ + - S X 
V333333 33 

4 ^ 2 2 52 2 4^ 2 32 ^2^ ^ 
+-X - - Yt:^ TT^ + -S^TT^ H - 

3 3 27 3 9 I u 

/16 ^2 32 ^ ^ 40 2 32 .\ t2 

+y^^-27^ 

\^ 9 9 27 9 ; s2 

+ f-X3-25X2 + -X2-2X7r2 + -X+-5X--7r2 + -52 
V 9 3 9 27 9 

+^ + 45] | + (3.34) 



1^ l^-X + r^ |^3X2-4Xy-14X + 3y^-6y + 2 7r2 + 4j ^— 

+ (ax^--x^y--x^ + -xy^ + -xy-ax^ + -x--y^-2aA - 

^333 3 33 j u 

+ (^X^ + 2X^Y + 2QX^ + ^XY -^XTi^ + ^-X-^-^A\ + L^t\. 



(3.35) 
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4. Summary 

In this paper we presented analytic expressions for the O (a^) QCD corrections to 
the 2 — > 2 gluon-gluon scattering process due to the interference of the tree-level 
diagrams with the two-loop graphs in the MS scheme. Throughout we employed 
conventional dimensional regularisation. 

The renormalised matrix elements are infrared divergent and contain poles down 
to O (1/e^). The singularity structure of one- and two-loop diagrams has been thor- 
oughly studied by Catani [21] who provided a procedure for predicting the infrared 
behaviour of renormalised amplitudes. The anticipated pole structure agrees exactly 
with that obtained by direct Feynman diagram evaluation. In fact Catani's method 
does not determine the 1/e poles exactly, but expects that the remaining unpre- 
dicted 1/e poles are non-logarithmic and proportional to constants (colour factors, 
TT^ and (s). We find that this is indeed the case, and the constant H^"^^ is given in 
Eq. (3.24). Its origin is in double emissions from the final state partons. It is related 
to that found for quark-gluon and quark- quark scattering in a straightforward way 
and therefore provides a very strong check on the reliability of our results. 

The pole structure of the two-loop contribution is described by Eq. (3.2) while an- 
alytic formulae for the finite part according to the colour decomposition of Eq. (3.29) 
are given in Eqs. (3.30) to (3.35). The one-loop contributions to the two-loop pole 
structure are expressed in terms of the one-loop bubble graph in D = 4 — 2e di- 
mensions and the one-loop box graph in D = 6 — 2e dimensions for which series 
expansions around e = are provided in Appendix A. 

The results presented here, together with those previously computed for quark- 
quark scattering [17, 18, 19] and quark-gluon scattering [20] form a complete set 
of two-loop hard scattering matrix elements for parton-parton scattering at O (af). 
They are vital ingredients for the ncxt-to-ncxt-to-lcading order predictions for jet 
cross sections in hadron-hadron collisions. However, they are insufficient to make 
physical predictions and much work remains to be done. A major task is to establish 
a systematic procedure for analytically cancelling the infrared divergences between 
the tree-level 2^4, the one-loop 2 — > 3 and the 2 — > 2 processes for semi-inclusive 
jet cross sections. Recent progress in determining the singular limits of tree-level 
matrix elements when two particles are unresolved [40, 41] and the soft and collinear 
limits of one-loop amplitudes [42, 35, 43], together with the analytic cancellation of 
the infrared singularities in the somewhat simpler case of e+e~ photon + jet at 
next-to-leading order [44], suggest that the technical problems will soon be solved 
for generic 2^2 scattering processes. 

A further complication is due to initial state radiation. Factorization of the 
collinear singularities from the incoming partons requires the evolution of the parton 
density functions to be known to an accuracy matching the hard scattering matrix 
element. This entails knowledge of the three-loop sphtting functions. At three- 
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loop order, the even moments of the sphtting functions are known for the flavour 
singlet and non-singlet structure functions F2 and Fl up to = 12 while the odd 
moments up to iV = 13 are known for F3 [45, 46]. The numerically small Np non- 
singlet contribution is also known [47] . Van Neerven and Vogt have provided accurate 
parameterisations of the splitting functions in x-space [48, 49] which are now starting 
to be implemented in the global analyses [50] . 

Finally, and most importantly for phenomenological applications, a numerical 
implementation of the various contributions must be developed. The next-to-leading 
order programs for three jet production that have already been written provide a first 
step in this direction [51, 52]. We arc confident that the problem of the numerical 
cancellation of residual infrared divergences will soon be addressed thereby enabling 
the construction of numerical programs to provide next-to-next-to-leading order QCD 
estimates of jet production in hadron collisions. 
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A. One-loop master integrals 

In this appendix, we list the expansions for the one-loop box integrals in D = 6 — 2e. 
We remain in the physical region s > 0, M,t < 0, and write coefficients in terms of 
logarithms and polylogarithms that are real in this domain. More precisely, we use 
the notation of Eqs. (3.27) and (3.28) to define the arguments of the logarithms and 
polylogarithms. The polylogarithms are defined as in Eq. (3.26). 
We find that the box integrals have the expansion 



Box^(K,i) = 



e^^T(l + e)r(l 
2sr (1 - 2e) (1 - 



+2e 



-2e' 




U^{x) - XU 



hu{x) + YU^{x)--X^U2{x) 



8 



-X^Y + -X^Y^ 
6 4 



16 



4 3 45 



+ (ii^i)|+e»(e3), (A.l) 



and 

Box^(s,i) 



e^^r(l + e)r(l-e)2 (n 
2uV{l-2e){l-2e) \li 



2\ <^ 



(X^ + 2mX) 



+e 



[-2U^{x) + 2XLi2(a;) - + 2YX^ - tt^X + 2(3) 



+i7r 2Li2(a;) +4FX-X2- 



7^ 



2Li4(;2) + 2Ui{y) - 2yLi3(a;) - 2XU^{y) + {2XY - X^ - 7r2)Li2(x) 

3 3 2 3 6 y 

+i7r(^-2Li3(x) - 2Li3(y) + 2yLi2(x) + ^X^ - 2X^Y + 3XY'^ 



-yr + 2C3 



(A.2) 



Box^(s,m) is obtained from Eq. (A.2) by exchanging u and t. 

Finally, the one-loop bubble integral in D = 4 — 2e dimensions is given by 



Bub(s) 



e^^r (1 + e) r (1 - e)' / /i 
r (2 - 2e) e \ 



2\ f 



(A.3) 
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